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ABSTRACT 



HST imaging and slitless spectroscopy are used to examine where the strong 
Lyman a emission escapes from the interstellar medium in the starburst galaxy 
Mkn 357. An Ha image shows that the ionized gas is mostly in a global 
wind, rather than associated with the individual star-forming regions seen in 
the optical and UV continuum. The Lyman a emission comes predominantly 
from the northwest side of the wind structure spatially, and shows a significant 
redshift relative to the optical lines. Both of these properties are signatures of 
seeing the line photons backscattered from the far side of a prolate or bipolar 
starburst wind, fitting both with escape calculations and evidence for winds in 
high-redshift galaxies with net Lyman a emission. Scattering is most important 
within this wind itself, rather than involving a surrounding neutral medium, 
as shown by the decreasing relative redshift of the line peak from 250 to « 30 
km s~^ between the center and edge of the detected emission. The Lyman 
OL emission exhibits strong asymmetry in comparison with both the starlight 
and Ha structures. These results add to the evidence that kinematics, rather 
than gas metallicity or dust content, are the dominant effect in determining 
which galaxies have strong Lyman a emission, and that powerful (and perhaps 
episodic) starbursts are common among Lyman-break galaxies as well as those 
discovered from Lyman a line emission. 



Subject headings: galaxies: starburst - galaxies: individual (Mkn 357) — 
ultraviolet: galaxies — galaxies: ISM 
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1. Introduction 

Lyman a emission plays an important role in the excitation of gas in H II regions. 
Paradoxically, this role is so important that very little Lyman a emission should escape to 
be observed. The mean free path of a Lyman a photon in a typical H II region is so small 
that many scatterings should occur before leaving the region (the standard "Case B" of 
photoionization calculations), with line photons escaping either when they are emitted or 
scattered near enough to the region's boundary, or shifted enough in frequency to increase 
the mean free path. This process, resonant trapping, means that Lyman a photons have a 
much greater probability of being absorbed by grains or converted to lower-energy line or 
two-photon emission than do photons in the adjacent continuum, so that we would expect 
to see virtually no line emission escaping from star-forming regions. Such calculations were 
set out in detail by Auer (1968) and Bonilha et al. (1979), showing both how effectively 
resonance-line photons are destroyed in static media, and how strongly expansion of the 
medium enhances their escape. Indeed, in the static case, Lyman a escape into surrounding 
neutral gas may serve only to backscatter the line photons into the ionized region again, 
further increasing the odds of their absorption by grains (Chen & Neufeld 1994). Prom 
these considerations, a simple model of gas in star-forming galaxies would predict emerging 
Lyman a only from a thin outer shell plus photons far enough into the line wings to reduce 
the incidence of scattering. 

Empirically, the situation in star-forming galaxies proves more complicated. Early 
expectations that protogalaxies could be found from intensely strong Lyman a emission 
(Sunyaev, Tinsley, & Meier 1978) were dashed first by the theoretical reahzation of 
how important scattering effects could be, and then by the widespread failure to find 
high-redshift line emitters as powerful as simple predictions based on star-formation rate 
and ionizing flux suggested. However, lUE spectra showed that some local star-forming 
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galaxies do exhibit significant Lyman a emission (Hartmann et al. 1988, Giavalisco et 
al. 1996, Keel 1998). The strength of the hne is not strongly coupled to metaUicity, so 
something more than destruction by dust grains in otherwise comparable systems must be 
taking place. HST spectroscopy of small regions within star-forming galaxies has shown a 
similar diversity, with some showing net emission equivalent widths as large as 70 A , and 
some showing no detectable emission superimposed on the stellar and interstellar-medium 
absorption (Lequeux et al. 1995, Thuan & Izotov 1997). It is especially noteworthy that the 
very low-metallicity galaxies I Zw 18 and SBS 0335-052 show only net Lyman a absorption 
(Kunth et al. 1994, Thuan & Izotov 1997), again showing that additional factors govern 
the observability of this line. 

Strong Lyman a emission is common from galaxies at high redshifts, even those with 
no sign of an active nucleus. In fact, Lyman a selection has proven to be an efficient way 
to find galaxies to high redshifts. This ranges from groups of objects at ^ = 2.4 — 4.1 
(Pascarelle et al. 1996, 1998; Francis et al. 1997, Venemans et al. 2002) to selection of 
galaxies at z = 4 — 6.5 (Rhoads & Malhotra 2001, Malhotra & Rhoads 2002, Rhoads et al. 
2003, Cowie & Hu 1998, Hu et al. 1999, 2003). Spectroscopy of large samples of galaxies 
selected from the Lyman break shows that about 1/3 have significant Lyman a emission 
(Shapley et al. 2003). Dynamical information has shown that some of the very large 
"Lyman a blobs" seen at z > 2 have properties consistent with very powerful starburst 
winds (Ohyama et al. 2003). 

Winds offer a promising environment for the escape of Lyman a, since the backscattered 
photons from the receding side of the wind have the maximal wavelength shift, and therefore 
the largest escape probability. This shift is enhanced when the wind scatters line emission 
from a central source such as a concentrated starburst or AGN, so that the mean wavelength 
of line photons is already redshifted before scattering in the wind. The blueshifted side 
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does not have this advantage, since the most favorable scattering geometry for escape 
points the photons away from the observer. Recent calculations by Ahn, Lee, & Lee 
(2003) underline this point, tracing the strongly preferred redshifting of observed photons; 
under ideal geometry, multiple line peaks corresponding to various numbers of scatterings 
might even be observable. Further work (Ahn 2004) shows that widely distributed dust 
may be required to fully suppress these multiple-scattering peaks, and yield the kinds of 
single-peaked P Cygni profiles seen in many high-redshift systems. Escape of Lyman a 
from gas in winds fits with the wavelength asymmetry often seen in spectra of both nearby 
star-forming regions and high-redshift star-forming galaxies. As stressed by Mas-Hesse et 
al. (2003) from analysis of HST aperture spectra, the dynamics of surrounding neutral gas 
may be much more important in the escape of Lyman a than the galaxy's dust content, 
and this fact could shape our selection of high-redshift galaxies from line emission strongly 
in favor of objects with global winds. This recognition further complicates the situation 
in star-forming galaxies; global starburst winds give ample amounts of outfiowing gas, but 
these winds can be so hot as to have a very small and patchy neutral component. 

Understanding where and how Lyman a emission escapes from star-forming systems 
is thus important for our interpretation of the early history of galaxies, as well as for the 
physics of the interstellar medium. However, spatially resolved measurements of Lyman 
a from nearby galaxies present formidable challenges. Geocoronal Lyman a emission 
compromises detections at very low redshifts, while foreground absorption from the Milky 
Way's own H I can affect measurements of galaxies with redshifts up to several hundred 
km s~^ along many lines of sight. The lack of suitable narrowband filters has made 
spectroscopic approaches the most effective way to seek line emission. A handful of objects 
were observed using lUE, as noted above, and further spectroscopy of compact galaxies and 
individual star-forming regions using the instruments on board HST has given detections 
and line profiles in several instances (Kunth et al. 1998, Lequeux et al. 1995, Mas-Hesse et 
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al. 2003). In each case, the emission shows a P Cygni profile, even when stellar background 
light is not an issue, implicating kinematic processes in what we can see of Lyman a. 

A scanning observation of M33 by Voyager 2 (Keel 1998), with one dimension of spatial 
information, showed diffuse Lyman a emission with no peak at the location of the very 
luminous H II region NGC 604, leading to the conjecture that the line arises in the diffuse 
interstellar medium rather than from individual H II regions. More recently, Mas-Hesse 
et al. (2003) have used HST spectroscopy to trace the extent of P Cygni structures and 
wavclenght shifts in star-forming galaxies. The new UV imaging capabilities brought by 
HST's Advanced Camera for Surveys (ACS) have been exploited by Kunth et al. (2004) 
to attempt imaging of Lyman a emission in nearby star-forming systems, an approach 
complemented by the spectroscopic technique used here. This paper describes the use of 
slit less Lyman a spectroscopy, with acompanying images in the continuum and Ha, to 
examine the spatial and velocity structure of Lyman a emission from a luminous starburst 
galaxy. Strong radiative-transfer signatures are seen, and the line emission comes from the 
far side of a starburst wind, probably a global superwind. This fits with less direct evidence 
from high-redshift galaxies, and suggests that such winds were commonplace in the early 
Universe. 

2. Observations 

2.1. Mkn 357 as a Lyman a Emitter 

Mkn 357 (PC 0119-1-229) is a promising system in which to study the escape of Lyman 
a emission from star-forming galaxies. Its optical emission lines show the line ratios of 
a starburst, at redshift z — 0.053 and thus well beyond contamination of its emission 
by Lyman a absorption from gas in the Milky Way. Ground-based imagery shows a 
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close companion and extensive, asymmetric Ho; emission. An archival lUE observation 
(Keel 1998) shows very strong Lyman a emission, guaranteeing enough signal to trace 
its spatial as well as velocity structure using HST. As one of the UV-brightest starbursts 
at such redshifts, Mkn 357 has also been a high-priority FUSE target to seek emerging 
Lyman-continuum radiation. 

I present here a four-way test for the roles of resonant trapping and velocity structure 
in the escape of Lyman a. A slitless spectrum around Lyman a provides both the spatial 
and spectral structure of escaping line photons. An Ha image shows the distribution of 
ionized gas, to distinguish the roles of radiative transfer and structure in the gas. The 
impact of dust extinction can be assessed from comparison of broad-band UV and red-light 
images. These continuum images also show where potential ionizing clusters are, and 
thus distinguish roles for an outflow wind and discrete H II regions in the ionized gas. A 
ground-based spectrum covering Ha sets the redshift zero point in assessing the velocity 
dependence of the escape of Lyman a. 

Using Ho=70 km s~^ Mpc~^ and a flat Universe with Qm = 0.3, Mkn 357 is seen at a 
luminosity distance of 237 Mpc (distance modulus 36.87), and the projected linear scale is 
1.03 kpc arcscc^^. This gives linear pixel scales of about 50 pc for the WFPC2 images (on 
the PC CCD) and half that for the STIS ultraviolet data. With this distance, Mkn 357 is 
fairly luminous for such a compact system, with Mb — —21.43, My — —21.67. This makes 
it about twice as luminous as the brightest system observed by Mas-Hesse et al. (2003), 
namely IRAS 0833-1-6517 at Mb = —20.8, which shows closely comparable metallicity. 
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2.2. Lyman a spectroscopy 

Lyman a was observed using the "slitless" spectroscopic mode of STIS on board HST. 
In practice, the wide 2" sht (aperture 52X2) was in place, to exclude geocoronal Lyman a 
and O I A 1304 emission while including virtually all of the Lyman a emission from Mkn 
357, which is observed at 1280 A . The telescope orientation put the sht in position angle 
—42°, roughly along the major axis of Mkn 357 and crossing some of the companion as 
well. The FUV MAMA detector and G140M grating gave a spectral scale of 0.0529 A 
pixel"^, with the spectral resolution set by the object's structure parallel to the dispersion. 
The sampling in the purely spatial direction was 0.029" /pixel. The spectral range was 
1245-1299 A . The exposure was 1560 seconds, cut short of the planned 4133 seconds by an 
electronics reset. The total Lyman a flux from the archival lUE observations (SWP 18199 
and 45198) was used as a consistency check on the STIS flux scale given this timing issue, 
using the pipeline conversion factor dijf2pt — 0.084 to derive a total in-sht flux from the 
values calibrated for surface brightness. The STIS total for a 200-pixel swath along the slit 
is 9.4 X 10~^^ erg cm~^ s~^, somewhat less than the lUE value of 1.6 ± 0.3 x 10~^^, but not 
by a large enough factor to require a calibration difference, since some additional Lyman a 
might faU beyond the 2" STIS wide slit. 

2.3. HST Imaging 

A far-UV image was also obtained with the STIS FUV MAMA, to complement the 
Lyman a data by tracing the morphology seen in hot stars, and particularly the presence 
and locations of massive ionizing clusters. The SrF2 fllter gave a passband about 120 A 
wide centered near 1460 A , with no signiflcant red leak. The total exposure was 2160 
seconds, with spatial sampling of 0.0247 arcseconds pixel"^. Both the STIS and WFPC2 
observations for this program were obtained on 12 September 2000. 
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The WFPC2 linear ramp filter (LRF) FR680N was used to isolate Ha emission, with a 
passband spanning 89 A FWHM. The central wavelength was set to 6922 A , 10 A redward 
of the Ha peak wavelength found spectroscopically. This offset, well within the FWHM of 
the passband and retaining about 90% of the peak transmission, allowed Mkn 357 to be 
observed on the PC CCD, with nearly Nyquist sampling of the image at this wavelength. 
A total of 1577 seconds' exposure was obtained in two images (with the second terminated 
by a guiding "loss of lock" shortly before its scheduled end), which were combined for 
cosmic-ray rejection. 

The continuum in red light was observed in a single WFPC2 exposure through the 
F702W filter, lasting 500 seconds. The pivot wavelength for this passband is 6917 A , 
close enough to redshifted Ha to make first-order color terms in the continuum subtraction 
irrelevant. Cosmic rays were removed by interactive editing, guided by the near-Nyquist 
samphng of the PC PSF at this wavelength, and the implausibihty of bright objects 
which do not appear at all in the narrowband image or UV STIS data. The galaxy was 
approximately centered in the PC CCD for this image, again taking advantage of its 
better sampling on a bright target. This pointing included the companion (or tail) of Mkn 
357, which fell off the chip in the Ha images. Alignment of the broad- and narrow-band 
images at the same telescope orientation should be a pure translation except for effects of 
optical distortion. Indeed, an image shift matches the structures in common as well as the 
difference image will show, without incurring the additional loss in resolution in resampling 
to correct for the small differential distortion between the two image locations. The 
continuum image was scaled for subtraction using the integrated Ha and [N II] equivalent 
widths from the INT spectrum (section 2.4), weighted by the LRF passband shape. To the 
approximation that all line emission within the F702W passband is included in the LRF 
filter, the two exposures contain enough information to separate the two contributions by 
linear combinations in which the constants are set by the relative filter widths and exposure 
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times. 

The resulting net Ha image shows ionized gas in a structure strikingly different from 
the starlight (section 2.2). This image was in turn subtracted from the scaled R image to 
give a nearly pure continuum distribution (except for the weaker [O I] and [S II] lines) for 
more accurate comparison with the UV starlight. 

The red-light image was aligned with the STIS UV continuum data using 9 compact 
sources (star clusters) which are bright in both passbands, with a fit accuracy per source 
of about 0.7 STIS pixel in each coordinate. This registration allows accurate ratio maps 
and photometry of the clusters, which were measured within circular apertures of radius 
0.1" using local background subtraction out to 0.25" radius. Their magnitudes were set to 
roughly the Johnson R zero point for the WFPC2 F702W data, and in the ultraviolet using 
the magnitude convention based on m = —2.5 log Fa — 21.1 where F^ is measured in ergs 
cm~^ s~-^ A 

2.4. Ground-based spectroscopy 

To define the systemic redshift accurately, and derive line ratios sensitive to the roles of 
stellar photionization and shock heating, I use a long-slit spectrum obtained with the 2.5-m 
Isaac Newton Telescope on La Palma, from December 1985. The 1.0" slit was oriented in 
position angle 135°, to include the companion object in the 1800-second exposure. The 
Intermediate-Dispersion Spectrograph was used with grating ruled at 831 lines mm~^, 
delivering a spatial sampling of 0.65" per pixel and spectral sampling of 1.1 A per pixel 
at 2.2 A FWHM resolution on a GEC CCD. The spectral range covered 6642-7289 A 
just including the shock-sensitive [O I] A6300 line observed at 6648 A . The systemic 
emission-line redshift from Ho: is z — 0.0531 ± 0.0001, with the uncertainty derived from 
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changes when different assumptions about the relative widths of Ha and the partially 
overlapping [N II] lines are used. The spectrum also shows a small net velocity gradient, 
of observed amplitude 53 km s~^ across the 18" span of detected emission, with the 
northwestern side receding. No line emission from the companion to the southeast was 
detected in this spectrum. 

This spectrum gives integrated equivalent widths (in the observed frame) of 242 A 
for Ha and 30 A for [N II]A6583. The profile of Ha is best represented as the sum of two 
components with Gaussian FWHM of 3.3 and 10.2 A (140 and 442 km s~^), with intensity 
ratio 1.9:1. In light of the WFPC2 Ha imagery, the broader component might reflect 
emission from escaping wind material. However, shock ionization in such an outflow is not 
important in the integrated hne ratios, with [O I]/Ha=0.02 and [S II]/Ha=0.16, neither 
of which is an unusual value for H II regions or young starbursts. This suggests that 
line emission from a wind is dominated either by in situ photoionization or line emission 
scattered from dust in the wind. 

The optical emission lines allow an estimate of the current gas-phase metallicity, using 
the traditional bright-line approach (Pagel et al. 1979, Edmunds & Pagel 1984, McCall 
et al. 1985). The INT spectrum was supplemented for this purpose by the blue portion 
of the Kitt Peak IIDS spectrum reproduced by Keel (1985). Using the relations from 
Pilyugin (2000, 2001) gives 12 + log (0/H) = 8.35, roughly half solar with the "new" solar 
cahbration. This is in the range seen nearby in such small starburst systems as NGC 1569 
(e.g., Devost, Roy, & Drissen 1997). 
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3. Starburst Properties 

3.1. Morphology and Interaction 

Ground-based imaging suggested that Mkn 357 was an inclined disk system with tidal 
tail and companion. These new data show a much more irregular structure, clearly not a 
normal spiral galaxy The companion, about 12" away in projection (12.3 kpc) does not 
show a central concentration, having a few diffuse bright knots within an amorphous glow 
that blends into the bridge connecting it to Mkn 357 itself. Two potential star clusters 
appear within the tidal debris around this bridge. These features are shown in Fig. 1, 
comparing UV and red images. Mkn 357 is not a very large galaxy; half the i?-band light 
is emitted within a projected radius of 0.9" (0.94 kpc), while the WFPC2 R image detects 
continuum light out to only a 6.5-kpc radius along the most extensive directions. Mkn 357 
falls in the range of size and luminosity often associated with blue compact dwarf galaxies. 

Within Mkn 357, both the UV and R images reveal a rich population of luminous 
star clusters, concentrated to the center, further evidence of its starburst nature. Local 
starburst systems are often taken as analogs of high-redshift galaxies. This is probably 
useful as regards high levels of energy input from star-forming regions, but will not generally 
be accurate in detail, since even rather quiescent galaxies today incorporate the heritage 
of a Hubble time's worth of gradual ISM enrichment and, quite likely, episodic infall of 
less-processed material. 

Further detail in the structure and history of Mkn 357 may be gleaned from a simple 
tomographic approach, suggested by the work of Katz-Stone & Rudnick (1997) on radio 
maps and foreshadowed by Zwicky's (1955) photographic decomposition of galaxy images. 
In a simple system consisting of two populations which differ in both color and spatial 
structure, an empirical separation may be performed by subtracting variously scaled 
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versions of an image in one band from an image at a very different wavelength, with the 
purest view of one population's structure found at the largest scaling factor for which 
significant oversubtraction does not occur. In realistic situations, a single scaling factor may 
be inadequate for a whole object, but useful insights can be gained from even approximate 
decomposition. For Mkn 357, the UV and R images decompose largely into two components 
differing in both structure and color (Fig. 2). 

The bluer component includes many of the prominent star clusters (section 3.3), and 
shows a rounder outer shape, with typical axial ratio 0.9. Existing ultraviolet observations 
of prominent starbursts (Marcum et al. 2001, Hoopes et al. 2004) suggest that this 
component may include scattered light from grains, perhaps including grains mixed 
with a starburst wind. The strong forward-scattering phase function of typical grains in 
the UV will weight this scattering toward those in the foreground of the stars (Witt & 
Gordon 2000), so this distribution is not necessarily representative of the overall grain 
distribution. The redder component is much more elongated, at least in projection, with 
typical axial ratio 0.3. Notably, there are several star clusters which appear only in the 
redder population, particularly at the southeastern end. The brightest of these might 
represent the nucleus of a pre-existing system which has seen strong tidal disruption, or be 
a remnant of a previous burst of star formation. This cluster is in fact redder than the rest 
of this component, leaving a hole in the blue component in a decomposition which works 
well almost everywhere else. 

3.2. Ha Morphology: a Stairburst Wind 

The continuum-subtracted Hq;-|-[N II] image (Fig. 3) shows structure quite different 
from that seen in the continuum. The dominant feature is an irregular spray of emission, 
less flattened in projection than the continuum light but more elongated than the blue 
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tomographic component, centered near the brightest star clusters. The brightest emission 
traces a U shape, with peaks close to (but not centered on) several bright clusters. This 
structure resembles the Ha appearance seen on a slightly smaller scale in NGC 1569 
(Devost, Roy, & Drissen 1997; Martin 1998). However, there is not a detailed match of 
continuum and emission-line features; in particular, the image does not show distinct H 
11 regions associated with most of the UV-bright clusters. This is a common situation in 
galaxies with strong and ongoing star formation. Numerous similar examples appear in 
the collection by Marlowe et al. (1997) and the extensive atlas by Gil de Paz, Madore, & 
Pevunova (2003). 

The prominent dark lane in Ho; does not seem to be associated with extinction, as 
traced either by a lack of star-forming knots or the UV/R broadband color. This again 
is reminiscent of NGC 1569, where a similar open structure appears in both optical line 
emission and the radio continuum, so it must reflect a genuine lack of ionized gas. In that 
system, this can be attributed to wind sweeping (Martin, Kobulnicky, & Heckman 2002), a 
picture which also makes sense in Mkn 357 given the other evidence for significant outflow. 

These characteristics suggest that much of the optical line emission arises in a global 
starburst (super)wind rather than from individual H II regions. This wind may be reflected 
in the broad (400 km s~^ FWHM) component seen in the integrated Ha profiles (section 
2.4). Such an outflow may be sampled (in part) by the velocity gradient in Ho;, which 
may have a counterpart component in Lyman a as well (section 4.2). In some starbursts 
with powerful outflows, the bipolar structure is evident only with detailed kinematics (e.g., 
Heckman et al. 1995), and in Mkn 357 the combined Lyman a asymmetries in intensity 
and velocity may offer the best evidence for a bipolar outflow geometry. Comparison with 
kinematic and morphological features seen in the gas of star-forming galaxies on a wide 
range of scales (Legrand et al. 1997, Gil de Paz et al. 1999, Silich et al. 2002, Lipari 
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et al. 2004) points up an ambiguity in the scale of outflows - is there really a single 
global superwind, or an amalgamation of expanding, overlapping shells around individual 
star-forming regions? These shells or bubbles would not necessarily escape the galaxy 
potential (as a wind, properly called, would), and would not individually extent over a large 
part of the galaxy image. Definitive discrimination between local and global outflows would 
require high-resolution kinematic data or X-ray observations. Weak evidence suggesting a 
global pattern may come from the degree or redshifting of the peak Lyman a emission and 
its systematic change with distance from the galaxy center. 

3.3. Clusters in the Steirburst 

Twenty bright clusters are well enough separated in these images for individual aperture 
photometry. These span a wide range in both brightness and color, with magnitudes 
R — 19.0 — 26.4, m(1400) = 19.2 — 26.0. The color range is not connected to magnitude, 
and spans m(1400) — R — —2.1 to -1-4.6. One luminous red cluster is much larger than 
the others and might represent the core of a former interacting component. Before any 
internal reddening correction, the brightest clusters reach Mr = — 18,Mi4oo = —17.6. This 
is already luminous enough to suggest that the clusters must be very young, since such 
systems as NGC 4038/9 have clusters only up to about absolute magnitude -14 in these 
bands (Whitmore et al. 1999). 

The combination of young clusters with a lack of discrete H II regions suggests that 
this system is in a stage comparable to the late-stage evolution of individual H 11 regions 
discussed by Bosch et al. (2002), when stellar mechanical-energy input dominates the gas 
dynamics. 
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4. Lyman a Emission 

4.1. Spatial distribution of Lyman a 

The distribution of Lyman a. emission is well traced in the direction along the STIS 
slit, and shows a striking asymmetry and structural differences from both Ho; and the 
starlight. Fig. 4 compares traces along the slit in these features and in the composite 
quantity of emission equivalent width. The emission-line flux was evaluated over a 2.1- A 
band, with continuum taken from a similar band just redward of the line location. This was 
intended to span the entire observed velocity range of line emission, but some variation in 
the absorption component may be included. 

The Lyman a emission distribution is very asymmetric, strongest on the northwestern 
side (which is receding according to the Ha and [N II] velocity gradient). The northwestern 
side of the emission encompasses 72% of the total, and emission can be detected three times 
as far from the core on this side (where the line trace disappears into the noise about 8" 
from the core). The asymmetry is similar to that observed in Haro 2 by Mas- Hesse et al. 
(2003). 

There is a narrow spike of line emission at the starlight core, but the rest of the Lyman 
a structure lacks continuum counterparts. Broadly, this means that it does not come from 
discrete H II regions around the UV-bright clusters, but from a more diffuse distribution 
of gas. Perhaps more interesting, structures traced by Lyman a and Ha emission do not 
coincide. Quite generally, this means that radiative transfer (likely modulated by the 
overall velocity field) dominates over either the distribution of ionizing sources or the gas 
distribution in determining where Lyman a escapes from the galaxy. This remains true 
even if nonradiative ionization (for example, by shocks within or at the interface of the 
wind) is important, since it is the difference in structure between these two recombination 
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lines that matters. 

The Lyman a/Ra flux ratio ranges from a low value of 0.56 (in a broad region 0.3-3" 
southeast of the core) to a high of 21, seen on the edge of the emission region 1.4-3.1" 
northwest. Redistribution of the line photons in wavelength during scattering is so effective 
that this line ratio can locally somewhat exceed its pure-recombination value. 

4.2. Lyman a velocity structure - scattering in a wind 

The Lyman a line shows prominent absorption as well as emission wherever the 
signal-to-noise ratio is sufficient, with a characteristic P Cygni profile. Much of the 
absorption must be due to radiative-transfer effects within the gas itself rather than to 
superposition on stellar photospheric absorption or absorption of background starlight. 
This is shown by the large velocity offset of the emission from the absorption core, the 
spatial persistence of the absorption well beyond the compact UV continuum distribution, 
and its change with distance from the center and thus with column density of scattering 
ions. The effects of this absorption on the line profile are difficult to extract without more 
complete information, since it may combine three effects - absorption of the continuum on 
stellar photospheres, absorption by foreground H I in the outskirts of Mkn 357, and an 
attenuation feature of continuum photons due to the large scattering cross section near the 
fine center (as described by Chen & Neufeld 1994). For this reason, I concentrate here on 
properties extracted directly from the spectrum. These are the results which are both most 
robust, and most relevant to the comparison with high-redshift galaxies. 

To extract peak wavelengths from the asymmetric emission line in such noisy data, 
a high-order polynomial was fit to the data points from 1275-1285 A ; orders from 21-45 
gave very similar results, with the values for 35 being quoted. The peak wavelength was 
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taken as the peak of this fit (noting that multiple peaks are seen at some slit locations). 
This smoothed fit incorporates all the statistically significant structure in the emission-line 
profile, as can be seen visually from examination of the residual pattern upon subtracting 
the fit from the original data (Fig. 5), which accords with the Poisson errors for the STIS 
data delivered by the calibration pipeline. 

The peak wavelength was measured by following the maximum value of the polynomial 
fits from one spatial increment to the next (using the IRAF routine reidentify). Errors 
for this peak estimate were assessed in a Monte Carlo manner, from 1000 realizations of 
appropriate Poisson errors added to the smooth fit for each line of the spectrum. The 
dominant velocity feature is a redshift of the peak velocity, greatest in the middle of the 
detected line emission rather than near the galaxy core. As shown in Fig. 6, the offset of 
the Lyman a peak from systemic velocity is as large as 250 km s~^ near the middle of the 
emission region, declining to 30 km s~^ toward the edges of the detectable emission. In 
these regions, the emission also shows evidence for multiple velocity components, so genuine 
kinematic structure may be playing a role in the inferred wind. As noted above, the Lyman 
a emission region appears very asymmetric with respect to both the galaxy's starlight 
and the recombination emission seen at Ha, which indicates that radiative-transfer effects 
operate in a very asymmetric way. As noted above, if most of the line photons originate in a 
small core region before scattering, the typical line Doppler shift along our line of sight will 
be twice the characteristic expansion velocity, due to the backscattering "moving mirror" 
geometry. 

This analysis has treated the STIS spectrum as a pure wavelength-sht position 
measure. The wide slit means that location of structures parallel to the dispersion may 
play a role. Such a role is unlikely to account for much of the structure in peak wavelength, 
simply because it would require a large U-shaped emission region with no counterpart at 
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other wavelengths. There is very hkely a spatial component to the apparent line widths, 
whose extent depends on the angular span of Lyman a emission. Much of the width of the 
Lyman a peak could arise in this way if the emission and UV continuum span comparable 
angles. However, the sense of the Lyman a velocity structure runs opposite to the direction 
that would result from resolved galaxy structure if the emission roughly follows the UV 
continuum light, suggesting that the velocity structure is unlikely to be compromised. One 
test of the extent of this smearing by spatial structure considers the widths of interstellar 
absorption features seen just blueward of Lyman a, of which the strongest unblended one 
in this spectral range is Si III A1206.5. At the nucleus, this line shows a FWHM of 2.5 A 
(or 590 km s~^), and is slightly deeper in core attenuation and stronger overall than seen 
in Haro 2 (Mas-Hesse et al. 2003). This line shows a core-to-edge change in radial velocity 
no greater than 60 km s~^ across 65 STIS pixels (1.7 arcseconds), a span hmited by the 
background starlight. Since the Lyman a centroid changes by nearly 300 km s"^ across the 
same locations, confusion by spatial structure of the galaxy is unlikely to be important. 

The sense of the velocity offset across the galaxy core is the same as that of the 
HcK velocity gradient, although the span and slope are much greater. The ground-based 
spectrum shows only a 50 km s~^ range in Ha velocity across the entire system. 

5. Discussion 

A comparison among Lyman a. Ha, and multicolor continuum data for the luminous 
and metal-rich starburst system Mkn 357 shows that the Lyman a emission is associated 
with a region which is very asymmetric about the galaxy center and likely associated with 
the receding half of a global wind. The remains some ambiguity about whether this might 
be overlappig winds from individual star-forming events, or a global wid which is no longer 
bound in the galaxy potential. Some evidence in favor of a global superwind comes from 
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fairly cohesive velocity structure in Lyman a emission across the galaxy, suggesting one 
mass of gas is responsible. 

These conclusions fit with the results of several previous studies (section 1) in showing 
that strong Lyman a emission is more a sign of kinematics (favoring outflow geometries) 
than metallicity (through the abundance of grains). Applied at high redshifts, this means 
that starburst outflows are very common among Lyman-break galaxies as well as the 
narrow- line Lyman a emitters which are common at z > 2. They also fit with the redshift 
behavior seen from various spectral features in Lyman-break galaxies (Franx et al. 1997, 
Adelberger et al. 2003 Steidel et al. 1998, 2003). As Steidel et al. (2003) note, in the 
presence of a strong outflow, interstellar absorption lines will be preferentially blueshifted 
compared to the systemic value, while backscattered Lyman a is preferentially redshifted. 
Since the strongest absorption features in the easily-sampled mid-UV range can have 
significant interstellar components, the systemic redshift can remain uncertain by several 
hundred km s~^. 

These results also reiterate that Lyman a cannot be used as a simple tracer of the 
star- formation rate, reinforcing the arguments of Kunth et al. (2004). They find knots 
and diffuse structures not obviously related to the SFR distribution, fitting in some cases 
with a wind picture. From their sample, ESO 338-IG04 in particular shows structures 
much like the wind in Mkn 357. Comparison with these more local objects points up the 
role of neutral outflowing gas, which may be only a minority constituent of hot, mostly 
ionized winds from objects undergoing powerful starbursts, and adds further complexity 
to the interpretation of Lyman a emission from high-redshift objects. Not only active star 
formation and outflows, but the reservoir of additional gas or multiphase nature of the 
outflow, may play significant roles. 
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Fig. 1. — Montage of ultraviolet and red-light images of Mkn 357 and its companion 
structure, rebinned to the same scale and orientation. The field show is 15.4 x 17.2 arcseconds; 
north is 6° clockwise from the top. A logarithmic intensity display was used to moderate the 
contrast; even so, the inner structure is best shown in the lower panels, enlarged by a factor 
4 and stretched for a wider intensity range. 

Fig. 2. — Tomographic decomposition of Mkn 357, from the UV and red-hght images after 
rebinning to the same scale and comparable resolution. The redder component is much more 
elongated than the UV-bright blue component of the galaxy, which may include scattered 
light from grains seen at a small angle in front of the bright clusters. For clarity, the intensity 
mapping in each case is logarithmic with a small offset from zero. The field is 6.0 x 4.9", 
and north is 6° clockwise from the top. 

Fig. 3. — The Hq;-|-[N II] structure of Mkn 357. Left, the net emission- line image 
after continuum subtraction based on the i?-band F702W image, displayed with an offset 
logarithmic scale. Right, the same image shown as contours overlaid on the red-light 
continuum image, corrected for emission-line contamination, to show the different structures 
in starlight and ionized gas. The continuum image has likewise been displayed with a 
logarithmic intensity scale. As in Fig. 2, the field is 6.0 x 4.9", and north is 6° clockwise 
from the top. 

Fig. 4. — Comparative distributions of Lyman a, Ha and the ultraviolet continuum along 
the slit. Upper left, trace of Lyman a emission intensity (in arbitrary units) as observed and 
with 0.5" boxcar smoothing. The line flux was summed over a 2.1- A band, with continuum 
subtracted based on an equal-width band just to the red. Upper right, the Lyman a slice 
is overlaid with the ultraviolet continuum, both from the spectrum near 1300 A and from 
the STIS direct image. The two contionuum traces were used to ahgn the data sets, and 
their agreement in structure confirms that the positioning of the STIS and WFPC2 data 
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is reliable. There is a Lyman a emission spike near the starlight core, but the Lyman a 
profile is much more extended and shows structure on both sides of tis core not related to 
UV starlight Lower left, the Lyman a. and Hct profiles are compared. Again, the detailed 
correspondence is poor, with peaks in the two lines not aligned. The increasing Lyman 
a/Ho; ratio with distance from the center is a sign of radiative-transfer effects. Lower right, 
the trace of Lyman a mission equivalent width, in the emitted frame, is shown with and 
without 0.5" smoothing. The increase in equivalent width going outwards indicates some 
combination of radiative transfer and ionization by non-local sources. 

Fig. 5. — Lyman a profile from the STIS wide-slit observation, displayed with linear intensity 
scales after 2 x 2-pixel block averaging. The data are shown as obtained and in the polynomial 
fit used to derive peak wavelength, along with the difference between them (amplified by a 
factor 10 in amplitude for display). The region shown spans 10.07 A on the horizontal axis 
and 7.5" in the spatial (vertical) direction, oriented so that the northwest end of the slit (in 
position angle —42°) is at the top. The absorption component of Lyman a is prominent left 
of center in each panel. 

Fig. 6. — Lyman a profiles along the slit. Each is centered at the indicated distance 
from the UV continuum peak, in arcseconds, with negative values lying to the southeast 
along position angle 138°, and represents data summed over the indicated slit length in 
arcseconds. The heavy line shows the polynomial fit used to derive peak wavelengths, and 
the difference between data and fit is shown at the bottom. For comparison, all spectra 
have been normalized to the same peak intensity. The inner regions show pronounced P 
Cygni profiles, with the peak velocity redshifted by fa 250 km s~^ in a velocity scale set 
from the integrated Ha and [N II] emission. Far from the core, the emission has multiple 
statistically-significant components, including both redshifted and near-systemic values; no 
blueshifted components appear. 
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Fig. 7. — Velocity peak of Lyman a, as defined by the fits in Fig. 5, traced along the slit 
track. Error analysis is complicated by the profile fitting, but the envelope of the velocity 
error may be seen from local excursions as the intensity drops. Adjacent points are averages 
of two original lines of the spectrum, and are almost completely independent statistically 
(except for a small crosstalk due to rebinning in accounting for distortion along the slit). 
The shaded curve at the bottom indicates the net intensity of Lyman a emission at each 
location, to show the alignment with emission structures and the relative importance of 
different velocity-offset regimes. The zero point was set by the mean redshift of Hq;-|-[N 
II] emission, and should be accurate to ±30 km s~^. The error envelope indicates ilu for 
the main peak, and was evaluated with a Monte Carlo approach based on the pixel-to-pixel 
noise, adding realizations of this noise to the smoothed fits and remeasuring the peak. The 
error depends on not only the line flux, but noth the width and profile of the emission. 
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